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I. INTRODUCTION 


This report contains the description of four computer codes developed for the HP-21MX 
minicomputer system housed in the Control System Branch of the Guidance and Control 
Directorate of the US Army Missile Research and Development Command (MIRADCOM). 
These codes (programs) were developed to aid in analyzing fire control systems being 
developed or proposed by the US Army. 


The names of the symbolic versions of the programs are SFRESP, $PLOT3, $PLOTX, and 
$DREAD. The names of the binary versions are FRESP, PLOT3, PLOTX and DREAD. 
Uses of these programs are described in the following sections, and listings of the symbolic 
versions are contained in the appendices. 


Both the symbolic and binary version of FRESP, PLOT3, and DREAD are stored on the 
flexible disc labeled SYSCPY. The symbolic version of PLOTX ($PLOTX) is stored on the 
M-II SYSTEM disc. It was stored here because the user must supply the subroutine DATA 
and then recompile and relocate the program. The libraries on the M-I] SYSTEM disc must be 
searched in the relocation process. 


All binary programs must be executed under the control of the M-II SYSTEM contained on 
the M-II SYSTEM disc. Even if relocated properly, none of the codes can be executed under 
an M-I SYSTEM since they were developed using FORTRAN IV. 


2. FREQUENCY RESPONSE PROGRAM 


The program FRESP was developed to aid in computing frequency responses of con- 
tinuous or sampled-data contro] systems. This program is stored on the disc labeled M-II and 
can be executed under control of the M-II System contained on the same disc. The symbolic 
version of this program is stored in the file $FRESP. 


The general continuous and sampled-data systems for which FRESP is applicable are, 
respectively, shown in Figures ] and 2. For the continuous case, the program computes 
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where s=jw.* There are |. forward paths between the input. R(s), and the output, C(s). Each 
forward path has I cascaded transier functions where I depends upon & , i.e., each forward path 
can have a different number of cascaded elements. Each transfer function has the general form 
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The user specifies d, L, the I's, n, m, the a’s and the b’s. 


For the sampled-data case, the program computes 
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s'= s+ jkw.. (4) 


In Equation (4), s = jw and we is the sampling frequency. Equation (3) is the truncated infinite 
series version of the sampled frequency response of the system given in Figure 2. In Equation 
(3) T is the sampling period (T = 277/ wg). All other variables are the same as in the continuous 


case. 


A. DATA INPUT FOR FRESP 


Assuming that FRESP has been loaded into the memory of the HP-21MX, it can be 
executed with the command 


ON, FRESP 


After issuing this command the input data that describe the system and control the flow of 
execution are requested by the computer. A summary of the required input data is given in 
Table / in the order that is requested. It should be noted that only numbers I and 15 are 
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column dependent. When requested for these, the user must start his response in column one 
to be assured of a proper interpretation, {t should also be noted that both of these have default 


values. 


All other inputs are free field. If more than one number is required, the numbers should be 
separated by commas or blanks. Each line of input is terminated with a carriage return €. 
For a free field input, if a slash (/) followed by a@Q is issued, the input can be continued to the 
next line. Ifa line of input has not been terminated, i.e., a (S) given, the input of the line can be 
restarted by issuing a RUBOUT. 


B. ILLUSTRATIVE USES OF FRESP 


In order to illustrate the pragmatic use of FRESP, two examples are presented in this 
section. The first example is a single-forward path, sampled-data system, and the second 
example is a two-forward path continuous system. 


(1) EXAMPLE I. Consideration is given to the system shown in Figure 3. The goal is to 
compute the sampled frequency response when T = 0.02513 for 0.1Sw<100. Ten terms are to 
be used in the series and 50 frequency response points are to be computed and plotted. The 
session with the computer for accomplishing this is shown in Table 2. 


On the first page of Table 2 the dialogue with the computer for describing the system and 
controlling the computational flow is given. After the data is printed, the computer queries as 
to whether or nota plot is desired. In this case a plot was desired; therefore, YES followed bya 
ERwas typed. (The YES is not shown in this case.) The screen was automatically erased, and 
the frequency response plot shown in Figure 5 was made. 


(2) EXAMPLE 2. In this example it is desired to compute and plot the frequency response 
of the system shown in Figure 43 for 0.01<w<99.9. The results of the session with the 
computer for this case are presented in Table 3 and Figure 6. 

C. LIMITATIONS OF FRESP 

The limitations of FRESP are: 


@ Maximum of 5 forward paths (channels) 


@ Maximum of 10 transfer functions per forward path 
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1+ 0.5s 


R(s) 


1 + .258 + 0.062582 


Figure 3. System for example 1. 


C(s) _ 


64 
64 + 3.28 + s2 


Figure 4. System for example 2. 
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@ Maximum of 30 transfer function numerator or denominator coefficients per forward path 
@ Maximum of 100 frequency response points. 


Any of these limitations can be changed by varying the appropriate dimensions in FRESP and 
the corresponding subroutines. 


3. 3-D PLOT PROGRAM 


The programs PLOT3 and PLOTX were developed for plotting a function of two variables 
in a plane. The function is mathematically assumed to be represented as 


z= z(x,y) (5) 


where z is a dependent variable and x and y are independent variables. The function z can be 
thought of as a three-dimensional (3-D) surface. 


To obtain a true perspective of this surface, it is necessary to project it onto a plane 
perpendicular to a line-of-sight (LOS). (See Figure 7). The LOS is uniquely defined by the two 
angles a and B. Angle B, called elevation, is between the LOS and its projection in the x-y 
plane, whereas a, called azimuth, is the angle between the negative y-axis and the projection of 
the LOS in the x-y plane. 


Although the two angles a and B uniquely define the LOS, they do not uniquely define the 
projection of the surface onto the plane perpendicular to the LOS. This is not done until an 
orthogonal coordinate system has been defined such that two of the axes lie in the plane. To 
accomplish this, the z-axis is projected onto the plane. This projection is then rotated in the 
plane by an angle y, called tilt. (Actually the surface will appear to be tilted by -y.) This forms 
an axis which is called W. The W-axis and the LOS are perpendicular and can be considered to 
form two axes of a right-handed coordinate system. The third axis of this system, called U, lies 
in a direction parallel to the vector formed by crossing a unit vector parallel to W and a unit 
vector parallel to the LOS. By projecting points from the surface upon the axes W and U, a 
true perspective is achieved when viewing from a LOS defined by a and B and witha tilt of y. 


Since there is no translation involved, projection is easily accomplished by representing 
each point on the surface as a vector. Then the dot products of this vector with unit vectors 
along the W and U axis give the respective coordinate points of the projected point relative to 
these axes. 


20 


ae ven 


-aj6ue yy pue ajBue SO? ey) Buyuyep 10) jeO}D}q “2 eaNBiy 


21 


nr renner ae oo 


4 
{ 
| 
‘ 
§ 
H 
3 
& M : 
; Yj 
f 
x ~ 
8 


A. USING PLOT3 AND PLOTX 


Assuming PLOT3 or PLOTX has been loaded into the HP-21MX, they can be run 


respectively with the commands 
ON,PLOT3 
and 
ON,PLOTX. 
The computer will respond by printing 
GIVE NO. Y’S, ALPHA, BETA, AND GAMMA 


The programs were developed for plotting families of constant y-curves, i.e., y is held constant 
while x is varied and z is recorded. Thus, the first number requested is the number of constant 
y-curves in the family. The other three numbers are the angles a, B, and y mentioned 
previously. These numbers are read using a free field format and should be typed on the same 
line with commas delimiting them and followed by a 


The programs access the data to be plotted through the user supplied subroutine DATA. 
This subroutine has four parameters in the call statement. The first is an integer variable that 
indicates the curve number, i.e., the first y, second y, etc. This number is supplied to DATA by 
PLOT3 and should not be changed in DATA. The second parameter is a 1024 element, real 
array. The x values for the curve being requested are to be stored in this array. The next 
parameter is a real scalar in which the value of y for the respective curve is to be stored. The 
third parameter is a 1024 element, real array for storing the z values of the curve being 
requested. The last parameter is an integer scalar for transmitting to PLOT3 or PLOTX the 
number of data points used to describe the corresponding curve. 


In essence the exchanges of data between PLOT3 or PLOTX and DATA are: 


e PLOT3 or PLOTX gives DATA the curve number 


22 


| 
| e DATA returns the x’s,ay value, the 7's and the number of data points, 1.e., the number 


of (7.x) pairs 
In developing the subrouting DATA the following rules should be adhered to: 


@ The x’s must be either in ascending or descending order, 1.e., random order of the x’s 


will cause incorrect plotting. 

@ The curves are numbered starting with the lowest value of y if -90°<a<90° and 
starting with the highest value of y if 90<a<270. (They are numbered starting with one 
and ending with the number of curves.) 

B. EXAMPLE USE OF PLOT3 OR PLOTX 


Consider the following function: 


s 
s+5 


ao he Rs 
ae = a + 36 


z,(s) = 0.5 (6) 


wheres is the complex variables = x + jy. Thus, the function z can be represented by Equation 
(5). It is desired to plot z for x on the range of -10<x<10 and y on the range -9.25<y<10.75. 
Twenty constant y-curves are to be plotted. For each y-curve [01 (z,x) pairs are to be 
generated. 


The SUBROUTINE DATA presented in the listings of both PLOT3 and PLOTX was 
written for generating the data to be plotted. This subroutine was attached to both PLOT3 and 
PLOTX; both were compiled and relocated. The session on the computer for loading and 
executing PLOTX is shown in Table 4. Initially, PLOT3 was occupying the computer 
memory. It was “OFF-ed” and PLOTX was loaded. PLOTX was then run with the ON, 
PLOTX command and it immediately requested the number of y’s, a, 8, and y. The line 20, 15, 
19, 0 was typed followed bya The execution began and finally the plot shown in Figure 8 
was produced. Using PLOT3 the same results were obtained. 


C. LIMITATIONS AND COMMENTS ON PLOT3 AND PLOTX 


The programs PLOT3 and PLOTX are, in essence, the same. The major difference is that 
PLOTX requires about 2048 less words of computer memory than does PLOT3. PLOTX was 
specifically developed to obtain this savings. 
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Both programs plot constant y-curves one ata time. For any y-curve the maximum number 
of (7.x) pairs must not exceed 1024 unless the dimensions on the z and x arrays are changed in 


the main program and in the associated subprograms. 
4. DIRECTORY READ AND MODIFICATION PROGRAM 


The program DREAD was developed for reading and modifying directories on the HP 
flexible discs. Such a program was needed to restore directories on discs that had been 
destroyed accidentally, e.g., changing discs without using the software dismount and mount 
routines. 


In order to use DREAD, it is loaded into the computer and run with the command 
*ON,DREAD 
The computer then prints the following message 
GIVE THE TRACK, SECTOR, AND NO. OF DIR. ENTRIES. 


At this point the disc has been dismounted (from a software point of view). The disc can be 
removed from the drive and a different one inserted. Then the user must type the track and 
sector numbers of the part of the directory to be read and the number of directories in the 
sector to be viewed and possibly modified. These numbers should all be typed on the same line 
and separated by commas. The directory on the HP discs is stored in TRACK 66, starting at 
sector 0. In TRACK 66, 239 directory entries can be stored. Each directory entry is 16 words 
long. Each sector on TRACK 66 contains room for 8 directories except sector zero which has 
room for only 7 directories (The first 16 words of sector zero are used to store disc format 
information.) The Kth directory entry will fall in the following sector 


Is = 14 *((K-1)/8) - 60 * ((K-1)/40) (7) 


where all computations are computer integer arithmetic and the files are numbered 
consecutively starting with two. This equation is only valid for TRACK 66. It would be 
unusual that the directory extended beyond TRACK 66. 


After the computer reads the track number, sector number and the number of directories to 
be displayed, it queries the user if he wants to make any changes. If he does, he can specify the 
word numbers to be changed, along with the changes. (For the format ofa file directory entry 
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see the RTE-M Programmers Reference Manual, Appendix F.) After a directory has been 
displayed and changes have been made, the changed directory entry is displayed. No changes 
have really been made at this point. When all the requested directory entries have been viewed, 


the process can be aborted by typing the number indicated. 
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The following is a file manager listing of the symbolic version of the program FRESP. The 
symbolic listing is contained on the SYSCPY disc in the file SFRESP. 
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The following is a file manager listing of PLOT3. It is contained on the disc SYSCPY and is 
located in the file $PLOT3. 
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“PRECEDING PAGE NOT FILMED 


The following is a file manager listing of the program PLOTX and the user supplied 
subroutine DATA. These are contained on the M-I1 SYSTEM disc, and both are located in 
the file SPLOTX. Ihe subroutine DATA generates the data for the plot shown in Figure 8. 
For plotting other data, the user must substitute DATA with the appropriate version and then 
compile and relocate, 
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The following is a file management listing of the symbolic version of the program DREAD 
This program is located on the dise SYSCPY. The symbolic version is in file $DREAD é 
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B Elevation angle of the LOS 
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